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Abstract 

Gohmann, Kliimper and Seel derived the multiple integral formula 
of the density matrix of the XX Z Heisenberg chain at finite tem- 
peratures. We have applied the high temperature expansion (HTE) 
method to isotropic case of their formula in a finite magnetic field and 
obtained coefficients for several short-range correlation functions. For 
example, we have succeeded to obtain the coefficients of the HTE of 
the 3rd neighbor correlation function < cr|a"| +3 > for zero magnetic 
field up to the order of 25. These results expand our previous re- 
sults on the emptiness formation probability [Z. Tsuboi, M.Shiroishi, 
J. Phys.A: Math. Gen. 38(2005) L363] to more general correlation 
functions. 

MSC: 82B23; 45G15; 82B20; 82B80 
PACS: 75.10.Jm, 02.30.1k, 05.70.-a, 05.30.-d 

Key words: correlation function; density matrix; high temperature expan- 
sion; nonlinear integral equation; quantum transfer matrix 
Report-no: OIQP-06-10 
to appear in Physica A 

1 Introduction 

Gohmann, Kliimper and Seel derived [I] (see also [21 El IH E]) a multiple 
integral formula of matrix elements of a density matrix of a finite segment 
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of arbitrary length m of the anti-ferromagnetic spin 1/2 XX Z Heisenberg 
infinite chain at finite temperature in a finite magnetic field by combining 
the quantum transfer matrix approach [6j- [TD] and the algebraic Bethe ansatz 
technique. Their formula generalizes the multiple integral formulae for zero 
temperature [TT1 [T2| [T3] to finite temperature case. This is a fundamental 
quantity since thermal average of any operators acting non-trivially on the 
segment of length m can be expressed in terms of their formula. Thus it is 
an important problem to perform this multiple integral and extract concrete 
numbers from it. Their formula contains an auxiliary function, which is a 
solution of a nonlinear integral equation. This nonlinear integral equation 
is essentially same as the one for the free energy [8], [9]. Thus to evaluate 
their formula consists of two non-trivial tasks: to solve the nonlinear integral 
equation and to integrate the multiple integrals. In our previous paper [Ti] , 
we applied the high temperature expansion (HTE) method to a multiple 
integral formula [3] of the emptiness formation probability P(m) for the 
XXX model, which is the probability of m adjacent spins being aligned 
upward, and succeeded to obtain the coefficients of -P(3) up to the order of 
42. As for zero magnetic field case, there is also numerical calculation [TS] 
for the multiple integral for m = 2, 3. The purpose of this paper is to expand 
our previous results on the HTE for P{m) [H] to more general correlation 
functions for the spin 1/2 isotropic Heisenberg chain in a magnetic field h. In 
section 2, we introduce the multiple integral formula of the matrix elements of 
the density matrix [T]. In section 3, we evaluate this multiple integral by the 
HTE method. Based on the result of the HTE of the density matrix, we will 
calculate the HTE of two point correlation functions (I3.1l) -( l3~3l) . In particular 
for zero magnetic field case, we have succeeded to obtain the coefficients of 
the HTE of a 3rd neighbor correlation function up to the order of 25 (cf. eq. 
(I3.10p ). Section 4 is devoted to concluding remarks. 



2 Integral representation of the density ma- 
trix 

The Hamiltonian of the spin- 1/2 isotropic Heisenberg chain in a magnetic 
field h is given as 

h = jjyi°* s +i + ^ + -~X>*, (2.i) 

3=1 3=1 

where crj, crj, cr| are the Pauli matrices which act non-trivially on the j-th 
lattice site in a chain of length L. Here the periodic boundary condition 



2 



&j + L — a j is assumed. 

Let us introduce 2x2 matrices: 

A = ( o o ) ' e i = (oo)' e 2 = (io)' e 2 = (oi) (2 " 2) 

These matrices are embedded into the space (C 2 ) lglL on which the Hamilto- 
nian ( 12.1 p acts: 



e 



_« = ® e ^® J® (i - J,) , (2.3) 

where 7 = + ef, a, /3 £ {1,2} and j G {1,2, ...,L}. The above Pauli 
matrices can be written in terms of these matrices: crj = &j\ + e^, crj = 

zejg — ieji, cr| = e ji — e i2- We a ^ so P u ^ °"/ = e ii anc ^ °j~ = e i2- 

Gohmann, Kliimper and Seel derived [1] an integral representation of the 

density matrix of the XXZ chain at finite temperature T. The isotropic 

(XXX) limit of their formula is given as follows. 

< e ift e 2/32 • • • e m/3m > - lim _ K 



j 

in 



x 

j=l«+l+i 



|a 



x det ^ A, rf r v( 2 - 4 ) 

where a(v) and G(v,£) are solutions of the following integral equations. 

i f\ h 2J f dy \og(l + a(y)) ( . 
loga(u) = -- + — — - / — — — - F , (2.5) 



T v(v + i)T 7r 1 + (v — y) 2 ' 



GM = _ 7 f i ; , + / * ' w g(».o , (2 . 6) 



(v - - £ - i) 7c 71 1 + (« - V? 1 + ' 

Here (tt„)™ =1 and (/3 n )J£=i are sequences of 1 or 2. We define the number of 
1 in (ce n )™-, as |a + | and the position n of j-th 1 in (a n )™ =1 as a\\ a + = 1, 

J 3 

1 < a± < < ■ ■ ■ < — m - We also define the number of 2 in 

(/?n)™=i as and the position n of j-th 2 in (/3„)™ =1 as (3~\ (3^- = 2, 
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1 < (3 X < j3 2 < • ■ ■ < P\p-\ < Tn. We shall put = . +1 for j G 

{1,2,..., \a+\} and $J = (3~_ {a+l for j G {|a + | + 1, |a + | +2, . . . , |a+| + 
The contour C surrounds the real axis anti-clockwise manner. a(v) is defined 
as a(f) = l/a(f). The emptiness formation probability is a special case of 
this density matrix element: P(m) =< ej\ej + i\ ■ ■ ■ ej +m ^i\ >. In this case, 
the multiple integral formula (I2.4p reduces to the one in [3J. 



3 High temperature expansion 

In this paper, we will calculate the HTE of the following two-point correlation 
functions for finite magnetic field h: 

< >=< ( e A ~ e jl)( e j+i{ + e i+i2)( e i+2i - e j+2 2 2 ) >, (3.1) 

< °]°U > 

=< (ejl - ejl){e j+1 \ + e j+l 2 2 )(e j+2 \ + e j+2 2 2 )(e j+3 { - e j+3 2 2 ) >,(3.2) 

< a t a j+2 >=< e i?( e i+il + ei+i2)e j+2 2 > • (3.3) 

At h = 0, one can express nearest neighbor and 2nd neighbor correlation 
functions in terms of P(2),P(3): < erjcrj +1 >= 4P(2) - 1, < crJerJ +2 >= 
8(P(3) — P(2) + |). Then we can immediately calculate the HTE of these 
correlation functions up to the order of 42 from the results in [T3]. We 
can also calculate the HTE of < a~^aJ +2 > up to the order of 42 from the 
relation < cr^o-J +2 >= ~ < a z ^o z ^ 2 > which holds for isotropic model (XXX- 
model) at h — 0. Moreover one can calculate the nearest neighbor correlation 
functions by taking the derivative of the free energy of the XXZ-chain with 
respect to the anisotropy parameter. On the other hand, < a z -a z - +z > can 
not be expressed only in terms of P(m) even at h = 0. 

At first, we will calculate the HTE of a{v) from the NLIE ([2J|. Note 
that this calculation is similar to the one for the HTE for the free energy [9] . 
We assume the following expansion for small \J/T\: 

J i ( d\{v) 2 | / \\ ,J \2 



a(v) =exp ( > ^a k (v) ( - ) ) =l + ai (v)-+ I + a 2 (v) j (- 




+ ( ^ + a x {v)a 2 {v) + a 3 (v)j (^) 3 (3.4) 

ai(-u) 4 ai(v) 2 a 2 (v) a 2 (v) 2 . N , N , A ,J, 4 

+ 2 + 2^ + ai ^) fl3 ^ + ^ + " ' • 

Substituting (13.41) into (I2.5p . and comparing coefficients of (J/T) m on both 
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sides, we obtain the following integral equation for each m (m G Z>i): 

M ^ h 2 \: f dy a -j^ + A m (y) 
a m {v) = — - + — — r d m ,i - / ~~ / 75-, (3.5) 



J f (f + i) / ' 7T 1 + (f — 

where A m (y) is made of {a k {y)}™~^: 

= 0, A 2 (y) = — — , A 3 (y) = , 

^4(2/) = ^92 '"" ^ - 1 

We can solve (13.51) recursively. The first few terms of {a m (v)} are 

h 2i h 2iv 

ai[v) = — — — — , a 2 {v) = — — — — + 



J v(l + v 2 Y zv ' J(l + v 2 ) (1 + v 2 ) 2 

a 3 (v) = 



aAv J 



J(l + v 2 Y 

2iv(3 + 3v 2 + 2w 4 ) - 3v 2 ) ih 2 v h 3 



3(1 + v 2 Y 3J(l+t» 2 ) 3 2J 2 (l + w 2 ) 2 12J 3 (l + w 2 )' 

(3.7) 



One sees that only a\{v) has a pole at the origin. One can also solve the 
integral equation ( 12. 6ft based on the results of (12.51) . We assume the following 
expansion for small \J/T\: 

00 / T\ k 

fl U • (3-8) 

fc=0 ^ ' 

Substituting (13.81) into (12.61) . we can obtain the coefficients. The first few 
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terms of {fl^f, £)} are 



9i{v,0 = t; ; — , ; „ 2w , ; — + 



i + v *){i + {v- o 2 )(i + e 2 ) 2 J (! + (« - o 2 ) ' 

i{2v-i)i 2 h(2 + 2v 2 - 2v£ + £ 2 ) 



^ (i + v 2 ){\ + { V - o a )(i + e) 2 2 j(i + ^xi + ( V - o z )(i + e) 

y 3 («, = -i( (6£ 2 + 8) t; 5 - f (9£ 2 + 10) t; 4 + 4 (2£ 2 + 3) t; 3 

(3£ 4 - 4£ 2 - 9) v 2 + 2 (4£ 4 + 9£ 2 + 6) v - £ (£ 4 + 3£ 2 + 3) ) 

/(3(^ 2 +i)>-o 2 +i) (e+i) 3 ) 

h({v 2 + 1) £ 4 - 2v {v 2 + 2) £ 3 + (2w 4 + 7w 2 + 1) £ 2 - 2v£ + 3v 2 - 1) 



2J(^ + l) 2 ((t;-0 2 + l)(e + l) 2 
h 2 (i£ - 2iv) h 3 



A J 2 (v 2 + 1) ((v - O 2 + 1) (£ 2 + 1) 24J 3 ((v - O 2 + 1) 
One see that only go(v,£) has a pole at u = £. Taking note on the above 
results and substituting (13.4j) and (I3.8P into f)2.4p . we can calculate the coef- 
ficients of the HTE for the density matrix in a finite magnetic field h only 
by taking residues at the origin. Then we can calculate the HTE of the two 
point functions (13. lft - (13.31) by using the result on the HTE of matrix elements 
of the density matrix. For example, the HTE of a 3rd neighbor correlation 
function at h = up to the order of 25 is given as follows: 



(3.9) 



< a*a* >= -t 3 - r + — z~ + 



58t 5 382t 6 545t 7 14473t 8 3227t £ 



j j+3 ' 15 15 63 315 405 

2715697t 10 3476257H 11 31700839* 12 1354469786* 13 
+ —— — + 



14175 155925 51975 868725 

405659182t 14 499339935195K 15 1774892665829* 16 

+ 363825 + 638512875 + 638512875 
337708981437007t 17 1373907665967541t 18 15824076710853959t 19 

10854718875 32564156625 + 168741538875 

2538703407655267579t 20 1840071379448207089t 21 
+- 



9280784638125 12993098493375 

945676556978993784629t 22 10941168707397851825299t 23 

- (3.10) 



714620417135625 16436269594119375 

251602687993890101833103t 24 4135206247498584771500543t 25 

+ — — + ■ 



49308808782358125 528308665525265625 
where we put t = J/T. We also calculated the HTE of < cr|cr| + 2 > f° r 
finite magnetic fields h = 2,4,6,8 up to the order of 30; < a^aj +3 > for 
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h = 2, 4, 6, 8 up to the order of 20; < o-+oj +2 > for h — 2, 4, 6, 8 up to the 
order of 30. In [TBI E] > t ne HTE of correlation functions for zero magnetic 
field were calculated up to the order of 19 by other method. It is remarkable 
that we can also treat nonzero magnetic field case and the orders of our HTE 
are higher than the ones in [161 E] • 

We have plotted the Pade approximations of the HTE of these correlation 
functions for finite magnetic fields h in Figures HE Once we introduce the 
HTE, we can also consider J < case by analytic continuation, though 
original multiple integral formula ( 12.41) was defined for J > case. We have 
also plotted quantum Monte Carlo simulation (QMC) data by Shiroishi [18] 
based on the open source software in the ALPS project [19], [20]. The system 
size of the simulations is L = 128. Our HTE results agree well with these 
QMC data except in the very low temperature regions. At least in the region 
where the Pade approximation seems to converge, < o"|cr| + 2 > for J > at 
h = and for J < 0, < <Tj<7j +3 > for J < and < a^aj +2 > at h = 
monotonously decrease with respect to temperature (here we omit a figure 
of < o"|o"j + 2 > for J < 0); < 0"J<r| + 3 > for J > at h = monotonously 
increases with respect to temperature. When h > in figures [T1[2lol peaks 
appear at non-trivial temperatures . 

4 Concluding remarks 

In this paper, we have evaluated the multiple integral formula of the density 
matrix by using the HTE method, and thereby obtained the HTE of several 
two point correlation functions. Together with our previous paper [T4] on 
P(m), we recognized that the HTE method is efficient to evaluate the mul- 
tiple integral formulae of correlation functions PQ-[3]. I n [23], the multiple 
integral formula (12.41) has been reduced to finite sums over products of single 
integrals for short segments of length 2 and 3. If a similar reduction is done 
for length 4, one will be able to calculate the HTE whose order is higher than 
the one in this paper. 

Acknowledgments We would like to thank Masahiro Shiroishi for discus- 
sions and providing us his unpublished QMC data [T8]. Mathematica Com- 
putation in this work was carried out at the Yukawa Institute Computer 
Facility. 
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Figure 1: Temperature T dependence of < cr|o"| +2 > for J > with a 
magnetic field h. We have plotted the Pade approximations of order [21,21] 
for h = and [15, 15] for h = 2, 4, 6, 8 ([n,<i] means that the numerator is a 
degree n polynomial of J/T and the denominator is a degree d polynomial 
of J/T). An exact value of < crjcr| +2 > at (T,h) = (0,0) was calculated in 
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Figure 2: Temperature T dependence of < 0*0*^ > for J > with a 
magnetic field h. We have plotted the Pade approximations of order [12, 13] 
for h = and [10, 10] for h = 2, 4, 6, 8. An exact value of < S/S^ >= | < 
cr|cr| +3 > at (T, = (0, 0) was calculated in [22] ■ 
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Figure 3: Temperature T dependence of < o"JaJ +3 > for J < with a 
magnetic field h. We have plotted the Pade approximations of order [12, 13] 
for h = and [10, 10] for h = 2, 4, 6, 8. 
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Figure 4: Temperature T dependence of < a^aJ +2 > for J > with a 
magnetic field h. We have plotted the Pade approximations of order [21,21] 
for h = and [15, 15] for h = 2, 4, 6, 8. An exact value of < c r |c r | + 2 >— 2 < 
2 > at (T, /i) = (0, 0) was calculated in 



a 3 a 3 
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Figure 5: Temperature T dependence of < (T^a^ +2 > for J < with a 
magnetic field h. We have plotted the Pade approximations of order [21,21] 
for h = 0, [15, 15] for h = 2, 4, 6 and [14, 16] for /i = 8. 
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